Observability of the total inflationary expansion 
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I consider the question of possible observability of the total number of e-folds accumulated during 
the epoch of inflation. The total number of observable e-folds has been previously constrained 
by the de Sitter entropy after inflation, assuming that the null energy condition (NEC) holds. 
The NEC is violated by upward fluctuations of the local Hubble rate H , which occur with high 
probability in the fluctuation-dominated regime of inflation. These fluctuations lead at late times 
to the formation of black holes and thus limit the observability of inflationary evolution. I compute 
the average number (AN) of e-folds accumulated during the last NEC-preserving fragment of the 
inflationary trajectory before reheating. This is the maximum number of inflationary e-folds that 
can be observed in principle through measurements of the CMB at arbitrarily late times (if the dark 
energy disappears). The calculation also provides a reasonably precise definition of the boundary 
of the fiuctuation-dominated regime, with an uncertainty of a few percent. In simple models of 
single-field infiation compatible with current CMB observations, I find (AN) of order 10^. This 
upper bound on the observable e-folds, although model-dependent, is much smaller than the de 
Sitter entropy after inflation. The method of calculation can be used in other models of single-field 
inflation. 



I. INTRODUCTION AND SUMMARY 

Inflation produces primordial metric fluctuations that 
may be observed indirectly through CMB measurements 
such as WMAP ^J. An observation of CMB at present 
corresponds to the measurement of the inflaton evolu- 
tion about 60 e-folds before reheating f3| . Assuming that 
CMB measurements will be possible at indefinitely late 
times, one might hope to deduce information about ar- 
bitrarily early stages of inflation. Of course, late-time 
acceleration (persistent "dark energy") can make it im- 
possible to observe CMB at very late times Q. There is, 
however, another limit on our ability to see towards the 
past. This limit is caused by violations of the null energy 
condition (NEC) during inflation. 

To make the following arguments more specific, let us 
consider a model of inflation driven by a canonical, min- 
imally coupled scalar field (f> such that the field evolves 
from a large initial value (/)in (perhaps near the Planck 
boundary (jjpi) to the reheating point (/) = 0*. A typical 
model of this type has the inflaton action 
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We assume that the inflaton potential V{cj)) grows mono- 
tonically with (f> and that the slow-roll approximation is 
vahd. In models of chaotic type, e.g. V{(j}) oc 0^", we 
then expect that ^ (p-m ^ ipPii ^tnd that </>in is deep in 
the fluctuation-dominated regime. 

The evolution of (p during inflation can be pictured 
as a random walk superimposed on a deterministic drift 
towards (j) = (j)^ [l-Q . The random walk can be modeled 
as "diffusion" in (p space with the diffusion coefficient 
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while the mean drift velocity is the time derivative 
the slow-roll evolution, 
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During the last stages of inflation before reheating, the 
trajectory (l)(t) is monotonic (0 < 0) with nearly unit 
probability, although there is always a small probability 
of an upward fluctuation ((/> > 0). On the other hand, 
in the fiuctuation-dominated regime an upward fluctu- 
ation of has around 50% probability because random 
fluctuations dominate over the slow-roll motion. 

We now note that an upward fluctuation of corre- 
sponds to an upward fluctuation of the local Hubble rate 
H and thus to a local violation of the null energy con- 
dition (NEC) 0, S]. A violation of the NEC due to an 
upward fluctuation of on a distance scale L ~ 
leads to the formation of a Hubble volume that looks 
like a black hole from the outside Q, and to a real 
black hole when the overdensity enters the local Hub- 
ble horizon at late times after inflation [l3, El. This 
can be understood qualitatively by noting that a Hubble- 
size region of approximately de Sitter spacetime with lo- 
cal Hubble parameter H has exactly the energy density 
A = -^MpiH^ that corresponds to a black hole with 
the Schwarzschild radius H^^. An upward fluctuation 
of H therefore leads to an increase of the energy density 
beyond the Schwarzschild limit. In this way, an NEC vi- 
olation in the far inflationary past will limit the lifetime 
of any future observers who may be trying to perform 
CMB observations at very late times. 

The main focus of this paper is an investigation of this 
limiting effect of NEC violations. For the sake of this con- 
sideration, I will assume that dark energy eventually de- 
cays, so that the late-time universe is not expanding with 
acceleration and the local Hubble radius grows without 
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limit, permitting (in principle) observations of the pri- 
mordial density fluctuations on arbitrarily large scales. 

It is interesting to determine the time range within 
which the NEC can be violated during inflation. For 
each random inflationary trajectory (j){t) there exists a 
well-defined time of the last NEC violation before reheat- 
ing, i.e. a time t-^ such that the NEC is violated around 
t = tN but then is not violated any more. The evolution 
of <j>{t) before < = iN is thus not observable even in prin- 
ciple. On the other hand, the evolution of (f) after t — tjsi 
is in principle observable: If primordial fluctuations on 
a distance scale L are produced at t > t^, these fluc- 
tuations will be observed through the CMB fluctuations 
at sufficiently late times when the scale L reenters the 
Hubble horizon (we are assuming that the dark energy 
does not prevent such observations). Therefore, it is only 
the statistics of the last NEC-preserving^ segment of the 
inflaton trajectory (j){t) that is — even in principle — 
accessible to observations. 

In the following sections we will compute (within an 
adequate approximation) the mean duration of the last 
NEC-preserving portion of the trajectory (j){t) until re- 
heating. When fluctuations are negligible, the field 
evolves according to the slow-roll equation and so 
the time needed for evolving from = </)i to reheating 
(taking into account that < (pi) is 
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(We write —v because the value of (j) decreases with 
time, so v{(j)) < 0.) This formula, however, cannot be 
used directly to compute the duration of the last NEC- 
preserving portion of the trajectory, for two reasons: 
First, the NEC-preserving portion of the trajectory de- 
pends on chance and is not confined within a fixed in- 
terval, say 0i]. Second, the initial stages of the tra- 
jectory (j)(t) belong to the fiuctuation-dominated regime 
where the evolution (/)(t) is not well described by the de- 
terministic slow-roll equation </) = v{(j)). 

Below we will compute the average duration (AiNEc) 
of the last NEC-preserving portion of the trajectory (j){t), 
where the average is performed over the ensemble of all 
comoving trajectories.^ Heuristically, we may attempt 
to determine a value cj) — <j)q such that the duration of 
the slow-roll trajectory between (j) = (j)q and </) = 0* is 
precisely equal to (AfNEc), i-e- we first compute (AiNEc) 



^ I talk about "NEC-preserving" rather than about "monotonically 
decreasing" trajectories because in models with several fields 
{4>l, ...,(p„), an NEC violation does not necessarily entail an up- 
ward fluctuation of a particular field 4>k{t). 

^ Thus we compute the "comoving" average rather than a "volume- 
weighted" average, which would require more complicated calcu- 
lations left for future work. See, e.g., Ref. [l2j for a review of 
comoving and volume-weighted averaging prescriptions. 



and then define </>„ such that 
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This value (pq can then be interpreted as the boundary 
between the fluctuation-dominated and the fluctuation- 
free regimes. One of the main results of this paper is a 
method of computing (A^nec); thus, the boundary (j>q 
between the fluctuation-dominated and the fluctuation- 
free regimes becomes a well-defined quantity. A precise 
definition of this boundary is relevant, e.g., for certain 
measure prescriptions for regulating eternal inflation 
IT^ as well as for attempts to count the observable degrees 
of freedom after inflation [l^, [l^ . 

The boundary between the fluctuation-dominated and 
the fiuctuation-free regimes can be characterized in terms 
of the dimensionless ratio of —v5t (the change of the field 
(j) due to slow roll during one Hubble timestep St = H^^) 
to \/2Ddt (the typical fiuctuation during the same time) , 
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is the first slow-roll parameter. We note that 6^(0) co- 
incides, for infiationary models of the type ([T|), with the 
inverse magnitude of the power spectrum of the primor- 
dial scalar fluctuation mode that crossed the Hubble scale 
at that time. 



Ps 



4^ 



V{4>) 



3£i M4, 62(0)- 



(8) 



(In the equation above, we neglected the slow-roll cor- 
rections since our final result will only depend logarith- 
mically on Ps-) The fiuctuation-free regime is character- 
ized by 6(0) ^ 1 and the fluctuation-dominated regime 
by 6(0) < 1 (i.e. by primordial fluctuations of order 1 or 
larger). However, this qualitative characterization can- 
not provide a sharply defined boundary value 0g sepa- 
rating the two regimes. 

We also note that the order parameter fi, which was 
used in Ref. [l3| to characterize the transition between 
the presence and the absence of eternal inflation, is re- 
lated to 6 by 
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The fluctuation-dominated regime was characterized by 
the condition 51 < 1 in Ref. [17[ i which again corresponds 
qualitatively to 6 < 1. 

It will be shown below that the rigorously computed 
value (AiNEc) can be approximated as 



(A^NEC) 
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where /(O; 4>) is the probability of the event that an in- 
flationary trajectory (j){t) starting at t = with the 
given value of will never violate the NEC. It will be 
found that / is close to being a step function, f{0;4>) ~ 
0{4)q — (j)), effectively cutting the integration at </> = </)g. 
The relevant value of (jfq will be determined from an ex- 
plicit analytic approximation for /(O;0). 

If we use the number of e-foldings, = Ina, as the 
time variable t, the same method yields the average 
number of e-foldings, (AiVNEc)i during the last NEC- 
preserving portion of the trajectory before reheating. Be- 
low we will compute (j)q and (AiVNEc) explicitly for infla- 
tionary models with a power-law potential. In a specific 
example, we will show that the average number of NEC- 
preserving e-foldings in the model with V{(j)) = ^^rn^cf)^ is 
of order 10^, if one assumes the model parameters that 
fit the WMAP data. In this model, the value of c/iq turns 
out to be such that b'^{(f>q) ~ 14. Thus, (f>q is well outside 
the fluctuation-dominated regime. 

These results can be compared with the^upper bound 
on the e-folds of inflation obtained in Ref. 15] . Assuming 
that the fluctuations are never dominant (equivalently, 
that the NEC always holds), it was found that the num- 
ber of observable e-folds of inflation must be smaller than 
the entropy Sds of the final de Sitter state after inflation. 
The latter is an extremely large number, of order 10^^° (if 
we use the current value of the dark energy density) . The 
present calculation shows that the limit on the number 
of observable e-folds is much more stringent. Therefore, 
the number of observable degrees of freedom, if it is ex- 
pressed through the entropy of the final de Sitter state, 
is in any case not directly related to the total number of 
observable e-folds of inflation. 

It must be noted that the value (AA^nec) is highly 
model-dependent. In the present paper, we perform com- 
putations only for models with V{4>) oc 0^" and derive 
a formula for (AA'^nec) [Eq. that shows a sensitive 
dependence on n. However, the method of calculation 
developed in this paper is sufHciently general so that the 
number of total observable e-folds can be computed in 
any other model of single-field slow-roll inflation. 



II. THE TIME OF THE LAST NEC VIOLATION 

The evolution of the inflaton 0(i) for a model of 
type ((TJ is a random process described by the Fokker- 
Planck (FP) equation (see e.g. [12,, JJ.]) 

<9tF(0,t)=L^P(0,O, (11) 
L^P = {d^DP) - vP) , (12) 

where the coefficients D{<j)) and v{(j)) were defined above. 
The FP equation is supplemented by appropriate initial 
and boundary conditions. The initial condition 

P(0,i = 0) =(5(0-0i„) (13) 



reflects the initial value of the inflaton fleld at i = 0, 
while the boundary conditions are imposed at the Planck 
boundary (f> — (p-pi (we impose the reflecting boundary 
condition) and at the reheating boundary (j) — (j)^: 

[d4,{DP)-vP]icj)put) = 0; d^{DP){cp,,t)={). (14) 

These equations describe the "comoving" evolution, 
i.e. P{4>, t)d(j) is the (infinitesimal) probability of having 
the value of the inflaton within the interval [0, -I- dcj)] at 
a fixed point in space and at time t. The "propagator" 
(i.e. Green's function) P{(j)i, (j)2,t) for the FP equation 
describes the probability density of reaching the value 
= 02 at time t starting with = 0i at time t — 0. The 
propagator is the solution of Ec^. (|lll) with respect to 02 
with the initial condition 

P(01,02,O)=5(0i-02) (15) 

and the same boundary conditions as the FP equation, 
namely Eq. (fT4|) . with respect to 02. The propagator is 
a function only of the time interval t since the evolution 
is invariant under time translations. The probabilistic 
interpretation of the propagator is that P(0i, 02, t)(i02 
gives the probability of the field value 0(t) being within 
the interval [02, 02 + ^02] at time t, while 0i = 0(0) and 
t are sharply fixed. 

A technical complication for the present considerations 
is that an equation describing the statistics of NEC- 
preserving evolution cannot be formulated as another FP 
equation, e.g. with modified coefficients. This is so be- 
cause the evolution strictly according to the FP equa- 
tion will violate the NEC at any time. Mathematically, 
FP equations describe a Brownian motion superimposed 
onto a deterministic motion, while the Brownian motion 
admits large velocity fluctuations at short time scales: A 
typical fluctuation (50 cx y/St over a time 6t produces a 

— 1/2 

velocity fluctuation cx {6t) , which is unbounded as 
St —J' 0. Thus, the description of the inflaton through 
Brownian motion effectively excludes the possibility that 
the trajectory 0(i) is strictly monotonic with < for 
any flnite duration of time. In reality, the mathemati- 
cal picture of Brownian motion does not hold for 0(t) at 
arbitrarily small time scales. The FP equations may be 
used to describe the evolution of only on time scales of 
order 6t ^ or larger. A formulation of the stochas- 
tic evolution restricted to the subset of NEC-preserving 
trajectories requires an averaging over such time scales. 

Therefore, one can formulate a statistical description of 
the subset of NEC-preserving trajectories only by using 
an equation that is nonlocal in time on time scales 5t, 
or nonlocal in on some relevant scale 50. Below we 
will derive one such equation and obtain its approximate 
solution. For now, we focus on determining the time 
of the last NEC violation, supposing that the statistical 
distribution of NEC-preserving trajectories is known. 

We consider the ensemble of comoving worldlines with 
inflaton trajectories <j){t) starting at t = with a flxed 
value 0(0) = 0in. The time is a random variable whose 
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distribution can be computed as follows. We ask for the 
probability Pr (^n < T) of the event < T, where T is 
a fixed parameter. The event < T means that the 
NEC holds for (j){t) after time T and until reheating but 
may be violated at any earlier time t < T. To compute 
Pr (tN < T), we split the random trajectory (j)(t) into two 
stages: The first stage is the evolution from = (f)in at 
t — Oto some intermediate value (/)t at time t — T\ during 
this first stage, the NEC may be violated. The second 
stage is an NEC-preserving evolution from </> = (pj- at 
time T until reheating at (/> = at some (random and 
not fixed) later time > T. It is clear that the last 
NEC violation happens before t — T ioi any trajectory 
consisting of these two stages, for any (jir and t*. On 
the other hand, trajectories with different values of (pT 
or are mutually exclusive random events. Therefore, 
we may simply integrate over all allowed values of 4>t and 

in order to compute the probability Pr [t^ < T). 

We will now compute the probability of the event that 
the trajectory (j)(t) has the two stages as just described. 
Since the first stage is not constrained with respect to 
the NEC, the evolution proceeds according to the FP 
equation. The probability of reaching an intermediate 
value (pT sX t = T IS thus given by the propagator 

P{(j)\n,<i>T,T)d(l)T. 

The evolution during the second stage needs to be 
NEC-preserving; in the single-field model we are consid- 
ering, this is synonymous with the trajectory (j){t) being 
monotonic. 

Since we are interested in the last NEC-preserving seg- 
ment of the trajectory (j){t) before reheating at = t/i*, 
we need to compute the probability density of reaching 
a fixed value (/> = (/>* at an unknown time , rather than 
of reaching an unknown value (j)2 at a fixed time t. Let 
us denote by P+{<j)o;t*)dt,, the probability of an NEC- 
preserving trajectory that starts at 0(i ^ 0) = po and 
reaches 4> = 4>^ within a time interval -f dt^]. (In 

this section, we will treat P and P+ as known; the neces- 
sary computations are postponed to the next sections.) 

Now we may express Pr [t^ < T) as the integral over 

and (f>T of the probability density 

Pr(<N < T:t^,,(l)T)dttd(j)T 

= P{(Pin,(l}T,T)d^TP+{(pT,t. - T)dU, (16) 

namely, 

Pr (^N <^)=J ^**J '^'t^T Pr (iN <T;U,(j)T) 

dU I d(f>T P(0i„, 0T, T)P+{^t; U - T). (17) 



Here and below, the omitted range of integration over (pT 
is from (/)pi to 

Once the probability Pr [t^ < T) is known, the proba- 
bility density p{t]<;) will be found from 



piT)^—PT{t^<T). 



(18) 



However, we will not proceed to compute p(iN) since our 
focus is on the duration At = t^, — t-^ of the last NEC- 
preserving segment of the trajectory (j){t). (As discussed 
above, the quantity At is observable in principle, while 
tN is not observable.) 

Let us denote by p{At) the probability density for At; 
by definition At > 0. It is more convenient to compute 
the generating function 



5(A) 



dTe^''p{T). 



(19) 



To assure the convergence of this integral, we will use 
g{X) only with A < 0. Once this function is known, we 
can compute the moments of the distribution p(At). For 
instance, the mean value of At and the dispersion a^t 
are given by 



(At 



NEC/ 
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(20) 



A=0 



al, ^ {Atl^c) - (A^NEc)' = 5"(0) - 5"(0). 

(21) 

The physical interpretation of (AtNEc) is the mean time 
spent in the last NEC-preserving segment of the trajec- 
tory 4>{t) before reheating, while dAt is the typical devi- 
ation from the mean among all the trajectories <^(t). 

In order to compute 5(A), we consider the joint proba- 
bility density of the time of the last NEC violation with 
the parameters t* and (pT', this probability density is 
found using Eq. p6|) as 



p{T;t,,ct>T) = ^Pr(^N < r;t„0T) 



?in, (PT 



,T)P+{^T,t,-T)]. (22) 



The value of g{X) equals the average of e'*'^* — e^^**^"^) 
among all trajectories that reheat at t = t* and contain 
the last NEC violation at t = T. Hence 
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ffT j dT dt^e^^*'-'^'>p{T;t^,(j)T) 
la Jt 

poo j*t^ 

i>T / dt* / dTx 
Jo Ja 

X e^^**"^^^ [P(0i„,</)T,T)P+((/.T,t* 



T)]. 
(23) 



Integrating by parts, we obtain 

*dre^(*--^'^ [P(0i„,</)T,T)P+((/.T,t* -T)] 

= P{(pin,C^T, t,)P+{(j)T,0) - e^'-PiPin, <pT , 0)P+ ((/)T , t* ) 

-f A / *dre^(*-~^)p(0i„,(^T,T)P+(0T,i* -T) 
Jo 

= Pi(pin,4>T,t*)S{(pT - (p*) 

+ X [ dTe^'^^'-^^P 
Jo 



(5(0T-0i„)P+(0T,t*) 
in, 0T,T)P+((/.T,i*-T). 
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Substituting this into Eq. (1^51) and simplifying, we find 



5(A) 
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+ X 



du[P{ 



DT I dT d{U - T) X 
/o Jo 



X e^(*--^'P(0i„,0T,r)P+(0T,i*-T) 

l-/(A;0in)+A / #T«'(0in,0T)/(A;</>T), 



where we defined the auxihary functions 

/(A;0)= / dre^^P+(0,T); 
Jo 



(24) 

(25) 
(26) 



So we will not actually need explicit expressions for the 
full distributions P(0in, T) and P+((^, r); it suffices to 
compute the functions / and 4'. 

As shown in Eq. (jAlip in Appendix El below . the func- 
tion 4'((/)in, (/)) can be approximated (up to slow-roll cor- 
rections) for (j) < by 



1 



(27) 



The dependence on the value of (/)in was omitted here 
because it is exponentially small as long as (f>in is within 
the diffusion-dominated regime. Thus we will omit the 
dependence on where appropriate. 

The function /(A; (j)) will be computed in Sec. HII] as 



/(A; 0t) = exp 



WiX;(l>)d^ 



(28) 



where the auxiliary function M^(A; ((>) is approximately 
found as the solution of Eq. below. We can then 
rewrite Eq. (1^^ as 



5(A) = 1 - /(A; 0i„) -f A 



d(t) 



/(A;0). (29) 



We note that /(O; 0t) is interpreted physically as the to- 
tal probability of never violating the NEC for a trajectory 
starting at <j) — (pT- If '/'in is in the diffusion-dominated 
regime, the probability /(O;0in) is exponentially small 
and can be neglected in Eq. (j^ . It follows that 
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NEC/ 
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/(O;0). (30) 



In the rest of the paper we will perform the calculations 
explicitly and show that the factor /(O; (j)) effectively cuts 
off the integration at a model-dependent value <f> — (pq^ 
which is in the regime where the diffusion is already small. 
A numerical calculation in a specific model of inflation is 
then given in Sec. IIVI 



III. DURATION OF NEC-PRESERVING 
TRAJECTORIES 

It is necessary for our purposes to compute the function 
in Eq. psp . which we denoted by /: 

/(A;0o)= / dTe^^P+(0o,T). (31) 

JO 

This is the generating function of the duration r of NEC- 
preserving trajectories starting at a given value (f> = (f>o 
at time t = and finishing at = 0* at an unknown 
time t = T. For instance, the mean duration of time un- 
til reheating among all the NEC-preserving trajectories 
starting at (j) — (f)Q is given by 



d_ 



In /(A; 00 ; 



(32) 



A=0 



As discussed above, we expect that the function 
f{X;(j)o) satisfies an equation nonlocal in (pQ. To derive 
this equation, we consider the change 6(1) = 4>{5t) — (piQ) 
of the (spatially coarse-grained) value of (p after a single 



Hubble time step 5t 
space. 



iJ at a given comoving point in 



,i) = v{cp)6t+iVWm. 



(33) 



where ^ is a normally distributed random variable. We 
denote for convenience by p(^) the probability density of 
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1 
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(34) 



The NEC-preserving property at the presently considered 
Hubble time step is equivalent to the condition 5(p < Q 



^2D{(P)5t 



2iJ2 • 



(35) 



(Note that the quantity h is always positive since H' = 
dH/d(p > due to the assumption dV/d(p > 0.) Since the 
probability P^(Pq,t) includes only trajectories that pre- 
serve NEC throughout their evolution, we must include 
only values of ^ such that ^ < b{<po) when we describe 
the Hubble time step leading from (po to (pQ + Sp. So we 
may express P+{(po, t) through P+icpo + 5(p, r — 5t) as 



P+{4>o\t) 



KM 



6t). (36) 



Here Sep = S(p{(po,^) under the integral is understood as 
a function of Using Eq. ((3T|) . we may now express the 
value /(A; (po) through the values of / at the next Hubble 
step as follows. We first integrate Eq. (1551) with e^'^dr 
from T = 5t to infinity and then exchange the order of 
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integrals and shift the integration variable r by St: 



St 
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^ / dT e^- / p(C) P+ {^Q + 5<j,,T- St) 

Jst J ~oo 

d^piO dTe^^^+'^^P+icj)o + Scj),T) 

-oo Jo 

'''^'^°^d^piOe^''fiX;c^o + Sc^). (37) 

-oo 



Note that the top line in Eq. ([371) is slightly different 
from the definition of f{X;(j)o): The integration proceeds 
from T — St rather than from t = in order to allow the 
subtraction t — St in the argument of P+ . The difference, 



St 



dTP+{<t>Q,T)^StP+{<t>Q,St), 



(38) 



is negligible as long as (t>Q is at least a few e-foldings 
away from reheating. This is so because P+{(f>o,St) is 
equal to the (exponentially small) probability of jumping 
from (f) = (jjQ directly to = 0* in one Hubble time St. 
Therefore, we may replace the top line in Eq. p7p by 
/(A; (j)o) and finally obtain the equation 



/(A;<^o) 



fc(0o) 



'0,0)- (39) 



This is the basic equation describing the function 
/(A; 0o); as expected, it is nonlocal in (j). 

It is not possible to approximate Eq. (p9|) by a diffu- 
sion equation (as is the normal procedure while deriving 
FP equations) because the integration in Eq. (p9| pro- 
ceeds over a (/)-dependent range. Rather than trying to 
solve Eq. (1391) directly, we will approximate the solution 
of Eq. by an adiabatic ansatz [Eq. pS)) below] . 

Up to now we have been using the proper time as the 
variable t. If a different time parameterization is desired, 
such as 



A{(p{t))dt, 



(40) 



where A((f)) is a known function, then the coefficients D, 
V, and St must be modified as follows. 



v=-, D^—, St^ASt, 
A A 



(41) 



while the dimensionless coefficient 6(0) is unchanged. For 
instance, passing to the e-folding time 



N 



Hdt 



(42) 



is implemented by choosing A{(j)) — H{(j)). Below we 
will compute (AA'nec) in a specific model of inflation by 
using this method. 



At the end of the calculation, we will only need to 
evaluate /(A = 0;(f>o)- As already mentioned above, 
/(O; 0o) is the fraction of trajectories that never violate 
the NEC among all trajectories (j){t) starting a.t (p = (j)o. 
We note that for 0o in the fluctuation-dominated regime, 
the probability /(O; 0o) rapidly decreases with growing 
00 because there is a signiflcant probability of violating 
the NEC at every Hubble time step at those 0o. On 
the other hand, the probability of violating the NEC in 
the no-diffusion regime is exponentially small, and hence 
/(O; 0o) is nearly constant and almost equal to 1 for 
00 in that regime. Therefore, we expect that /(A; 0o) 
has exponentially strong dependence on 0o- Moreover, 
P+(0», r) — 5{t); this can be shown by considering 



/•oo 

/ P+(0»,r)dr-l, 
Jo 



P. 



Oforr > 0, 



(43) 
(44) 



which holds because trajectories starting with = 
immediately reheat and have zero duration. Therefore 



/(A; 00 



/ e^^P+(0„r)dr = 1. (45) 
Jo 



Motivated by these considerations, we represent the ex- 
ponential behavior of /(A;0o) and the boundary condi- 
tion (|^5]) by the ansatz 



/(A;0o) = exp 



00 



W{X; (t))d(t) 



(46) 



where PF(A; 0) is a new unknown function such that 
1^(0; 0) > 0. We then divide Eq. ^ through by 
/(A; 0o) and expand /(A; + (50) to flrst order in (50: 



(^^^(06 exp 



4>o+S4> 



W{X; (j))d(j) 



<Po 



d^piOe'^'exp [-W^(A; 00)50(00,0] 

l-b 



_ ^{\~vW)St 



d^ 



' — oo V 
^[X-vW+DW^)St 



2ti 



exp 



-^-VWsiWS, 



1 1 Jb + WV2DSi\ 



(47) 



where we suppressed the argument 0o in the last line, b = 
6(0) was deflned in Eq. (1351) . while erfcc is the standard 
error function 



erf X = 



dt. 



(48) 



We now note that Eq. (j47|) does not contain derivatives 
of W; this means that we are using an adiabatic approx- 
imation where W^(A; 0) is assumed to vary slowly with 0. 
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Thus W{X; (f) is the unique real root of the transcenden- 
tal equation 



b + WV2DSt 

71 



exp [{vW - DW^ -X)St] + ierf 

(49) 

such that W(A; 0) > for A = 0. 

The coefficient &((/)) measures the influence of quantum 
fluctuations on the evolution (j){t). It is possible to obtain 
approximate solutions of Eq. in the cases & » 1 (a 
nearly fluctuation-free regime) and 6 <C 1 (a fluctuation- 
dominated regime). To simplify calculations, we pass to 
a new dimensionless variable r by rewriting Eq. (j49|) as 



r 



In 



ierf-^ 
2 V2. 



L{r) 



&2 

2 



X6t, 



(50) 



where r(A; 0) is related to W(A; (j)) by 
We solve Eq. ((50|) by using the inverse function ^^^(x), 



which then gives the solution W{X;(f>) through Eq. (|5T 
as 



iy(A;0) 



L-i(i62_A(5t)-5 



(53) 



Derivatives of W{X; </>) with respect to A can be ex- 
pressed through W{X;<j)) by computing the derivative of 



L'{r) = r + 



2 e-^^ 
.1 + erf^ 



L-\x) + 



(54) 



(55) 



and using Eq. ([55]) . For instance, we find 

dX V2m L'[L-HU^ - X6t)] 



2D e-i'' +^^' 
1 w + 



Further derivatives with respect to A can be obtained 
similarly. Since we will ultimately compute the generat- 
ing function .g(A) and its derivatives only at A = 0, it is 
sufficient to set A = in what follows. 

We also note that the function L{r) is "universal" in 
the sense that its definition does not depend on the infla- 
ton potential V^cf)). It is therefore useful to approximate 
the inverse function L~^{x) semi-numerically. We first 



consider the asymptotic behavior of L~^{x) for large x. 
Using the well-known asymptotic representation of the 
error function. 



erf X = 1 
we obtain 



1 



xwn 



exp I 



{l + 0{x~^)), x-^+oo, (57) 



and hence 



2a;V27f 



[l + 0(r-2)] , -f-cx), (58) 



[l + 0(x"i)] , .T^+cx). (59) 



This asymptotic formula allows us to obtain the approx- 
imate solution M^(0; 0) for the case 6 ^ 1 as 



(51) T4^(0;(/)) = 



i-i(i62)_5 l + 0(6-2) 



y/2DSt 



e 2 



/2^ b^V2DSt 



(60) 



The solution in the opposite regime 6^1 can be found 
by starting with the numerically obtained value 



L-\0) « 0.7286 EE ro 
and by expanding L^^{x) near a; = 0, 
d 



(61) 



L (x) ~ tq + X 



dx 



[L-']+Oix') 



ro + x ro 
Hence for 6 <C 1 we have 
To - 6 



1 



+ 0{x^), x^O. (62) 



2" 



0{b^ 



w{Q- (t>) = L , r"^ , b ^ 0. 

V2D5i 

(63) 

We have thus obtained the solution W{0; (j)) in the two 
opposite regimes. An approximation that holds uni- 
formly for all positive b can be obtained, if desired, by 
matching the asymptotic expressions near & = and 
b — 00, for instance, using the following interpolating 
function, 



(56) W{G-(j)) 



& + 0.7194 



^/2D5i 63^2^ + 1.80362 + 2.7286 + 0.9874 

(64) 

Numerical verification shows that this function approx- 
imates W^(0, (f)) to within about 2.5% relative precision 
for all & > 0. (We note that W is model-independent 
only as a function of b and \/2DSt, while b{(j)), D{(j)), and 
6t = H~^{(j)) of course depend on the chosen model of 
inflation.) 

However, it turns out that the approximation in 
Eq. (|60l) . which holds in the fluctuation- free regime, is 
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sufScient for our present purposes. Let us derive the cor- 
responding approximation for the function /(O;0o)j 



/(O; (/)o) = exp 



W{0;(l))dcl) 



(65) 



assuming that is such that 6^(</>o) ^ 1- Since W{0; cj)) 
is quickly growing with 0, the integral under the expo- 
nential above is dominated by the upper limit, so we can 
use the asymptotic estimate 



W{0;(f>)d(f>: 



4>o 



[i + o{b-')] 





db' 


-1 


H 53 







(66) 



In deriving this estimate, we neglected terms of order b^^ 
as well as derivatives of H and b, since these are merely 
slow-roll corrections. 

For (jjQ near reheating, we have /(O; 4>q) ~ 1 with expo- 
nential precision. The value of (jjq at which /(O; (j)q) first 
drops to exp(— 1) can then be found as the solution of 
the equation 



63 



db_ 
"d4> 



(67) 



0=05 



where we need to substitute bq = h{(j)q). This can be in- 
terpreted as a closed-form equation for bq if we express 
db/d(f) and H{ipq) as functions of bq. A numerical calcu- 
lation needs to be performed to solve this equation for bq 
in a particular inflationary model and to check that the 
resulting value of bq satisfies 6^ » 1, which is required for 



the validity of the approximation used to derive Eq. (1671) . 
(For instance, the calculations in the next section show 
that 6^ « 14 ^ 1 for the inflationary model with the 
potential V{<j)) oc (f>'^.) 

We have thus determined (j)q such that /(O; (pq) = . 
The function /(O; (/>) has a sharp dependence on and 
interpolates from 1 to within a narrow interval around 
(f) = (pq. To estimate the width of this interval, let us find 



the value (pq such that /(O; 
be determined as a solution of 



e ^. This value can 



27r 



Hi 



if') 



63 



db 



2. 



(68) 



Then the width of the interval in 6„ can be estimated 



— bq. Since the exponential above is the 



as 6bq = b'"q^ 

fastest-varying function of 6, to first approximation we 
have 



6^(0(2) 



)«62(,/,,)-21n2. 



(69) 



(Numerical calculations show that this is an overestimate 
of Sb by about 20%.) This leads to a change in the value 



of (pq that can be computed through 
dhib' 



6b 
~b 



d\n( 



Sb 
T 



-6^2 In 2. 



(70) 



Since 6^^ ^ 1 while the logarithmic derivative 
9 In 6/9 In (/i is not large, we find that 6(pq /pq ^ 1 . Hence, 
the function /(O; 0o) has the effect of a cutoff near p = pq 
when integrated with a slowly- varying function of (f) such 
as l/v{(p), as required for Eq. ([50]) . 

Below we will replace integrations with the factor 
/(O; (f)) by integrations with the upper limit p = pq. This 
approximation introduces a certain error; to estimate the 
effect of this error on the calculation of (AiVNEc), let us 
find the number of e-folds in the slow-roll trajectory be- 
(/'g'^^ and (p = (pq: 



tween 



6N„ 



^ H{P)dP _ Hi(Pq)S(Pq 

H6b 



-V{(p) -V{(pq) 



-bv 



din 6 



77 In 2 



dhib 



To estimate d\nb/d(p, we use Eqs. © and find 

dln6 _v' ID' IH' _ £20? 
d(P ^ ~^ 215 ^ 2IT ~ ^Mpiy/eT' 

where S2 is the second slow-roll parameter, 

■y/2 y/n 



£2 



V 



We also have 



-V 

H 



PI 



Therefore we obtain the estimate 

In 4 



(71) 



(72) 



(73) 



(74) 



(75) 



This estimate is important because it displays the er- 
ror inherent in the definition of the boundary of the 
fluctuation-dominated regime. Below we will check that 
this error is acceptable when determining the average 
number of observable e-folds. 

Let us summarize the calculations presented so far. We 
have derived an estimate of the mean time (A^nec) of 
the last NEC-preserving portion of the inflationary tra- 
jectory: 



(At 



NEC/ 



/(O;0) 



(76) 



where the value of (pq is determined from Eq. (|67|) . We 
have introduced a simple approximation where the func- 
tion f{Q\(t>) is replaced by a step cut-off aX p = pq] the 
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error of this approximation is expected to be small. It is 
possible, in principle, to determine the function /(O; (f>) 
numerically and thus to obtain a sharper estimate, as 
well as to compute the standard deviation cfan of the 
number of observable e- folds using Eq. (|2ip . However, 
we expect that the standard deviation (tan will not be 
larger than the width of the function /(O; (f>) around the 
point (f) = <j)q, which is of order 5(j)q as estimated above. 
Therefore, it will be sufficient for the present purposes 
to use the estimated width of the function /(O; (j)) as the 
statistical uncertainty in (AiNEc)- 



The number A''e of inflationary e-foldings accumulated 
between some value (j) = <f>i until reheating is estimated 
(assuming (/>! as 



Ne{4>l) 



H{(j))d(j) 2TT(j)\ 



-v{4>) 



(84) 



The squared amplitude of scalar primordial perturbations 
generated at (/> = is given by WMAP observations as 
Ps w 2.3 • 10"^ (we use the data from Ref. '^]). As- 
suming that this amplitude is generated at e-foldings 
before reheating (below we will set « 60), we find 



IV. EXAMPLE: INFLATION WITH A 
POWER-LAW POTENTIAL 



PI 



2n 



> 



(85) 



We now perform specific calculations of 4>q and the av- 
erage number of NEC-preserving e-folds, (AA^nec)! for a 
model of single-field inflation of type ([ij with the poten- 
tial 



1/(0) 



A/, 



PI 



(77) 



This model can fit the current observations when n = 1 
or n = 2 (see, e.g., [l9|). We use the formalism developed 
in the previous sections for computing (AfNEc), except 
that we divide St, v{(j)), and D{(j)) in every formula by 
the factor -ff (0) in order to pass from the proper time t 
to the e-folding time N . 

For this model, we find in the slow-roll approximation 




n+l 



(78) 
(79) 
(80) 



The slow-roll parameters (computed as functions of 
through the potential V) are 



M2 T/'2 

£1 = 



£2 = 



167r y2 

ya yll 



477 



V 



n M^j 



(81) 
(82) 



Reheating is assumed to happen at (f> — (f>^ with £1(0*) 
1, which gives^ 



and then, using Eq. ([5]) with 



A = 



327r 



n^Ps 



2n 



)i, we get 

n+l 

< 1. 



(86) 



Substituting this value of A and Eqs. (I7Hl)-(|Hni) into 
Eq. (j67p . we can obtain an implicit equation for the value 
(j)q. However, it is more convenient to express (pq through 
hq = h{4)q) using Eqs. (|80l) and 



PI 



p2(n + l) 



(87) 



and to derive a closed- form equation for h^, 



" n+l Ps 



-21nAr,-21n 



■K {n+l) bq 



3n + 5 

+1 



This equation is in a form that can be solved numerically 
by direct iteration. We use the values Ps = 2.3 • 10~^ 
and Ne = 60, while n can be 1 or 2 [T^, and obtain 



13.9, 
8.9, 



2. 



A rough order-of-magnitude analytic expression for hq 
can be obtained by using only the first term in Eq. 



+ 1 Ps 



40 



1 



(89) 



However, this overestimates 6^ by about 50%. 

Let us now compute the average number of e-foldings 
during the last NEC-preserving part of the trajectory, for 
potentials F(0) oc . According to the results of the 
previous section, we need to integrate the e-foldings until 
the value 0g. Using Eq. ([87)) . we find 



PI 



47r 



(83) 



(A7V] 



NEC/ 



4,, ~v[(t)) 



2-K 



n Ml^ 



^ The value of is only an estimate because it is computed in 
the slow-roll approximation, which does not hold near reheating. 
However, our results are not sensitive to the precise value of 0* . 



= 0[l)Ps 



21nPe 



1 



(90) 
(91) 
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In the last line, we substituted for bq the rough esti- 
mate (j89p merely in order to obtain a simpler analytic 
expression for (AA'nec) a-s an explicit function of the 
parameters. We use the more precise Eq. ([TO]) for com- 
puting the numerical answers. With the value n = 1, 
which is preferred by observations, wc find 



NEC/ 



3.4- 10^ 



(92) 



With n = 2, we get (AiVNEc) ~ 2.2 • 10*. 

A curious coincidence is that Eq. (|90| can be expressed 
through the slow-roll parameter £2 very simply as 



(AiV] 



1 



NEC/ 



£2 



Here we need to compute £2(0) at the value (j)q 
mined through bq as b{(f>q) = bq. 



(93) 



deter- 



The value of (pq that corresponds to the obtained value 



of bq can be expressed as 



^ n N 2- 



PI 



2tt 



= — (ATVnec) 



(94) 



1 we have 

13] 



5.4-10'^Mpi, and for n = 2 



Thus, for n 

we have 4>q « 7.0 • lO^'Mpi. 

We also need to check whether the intrinsic error of 
the present approximation, as given by Eq. ([75]) , is small 
in comparison with the mean value (AiVNEc)- We get 



ln4 

SNq = , , 



In 4 



(ATVnec) ■ 



(95) 



With the numerical values used above, we find SNq ~ 10 
for n = 1 and SNq ~ 10^ for n — 2. The relative error of 
the approximation is given by ln4/(7r&g) and is about 3% 
for n = 1 and about 5% for n — 2, which is acceptable 
for the purpose of our estimates. 
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Appendix A: Solving the stationary FP equation 

The time integral of the propagator of the FP equation, 



?in, (PT 



dTPi 



?in, VT 



,T), (Al) 



can be computed in closed form in one-field models [2C 
One integrates Eq. (fTTj) in time and finds 

P{(f>o,cb,t)dt^ dtP{(f>o,(f>,t)dt 



P(0,oo)-P(0,O) = -J(0-0o) 

(A2) 



since P{(j>, 00) = 0. Hence, the function ^{(po, (f>) satisfies 
the equation 



(A3) 



with the same boundary conditions in as the distribu- 
tion P{(po, (j),t). In other words, ^ is the Green's func- 
tion of the stationary FP equation. The function can 
be computed by integrating Eq. (IA3|) . First, one obtains 



84, (D*) - V* 



bo) -t- Ci = 0, (A4) 



where the integration constant Ci is expressed using the 
boundary condition (|14|) at = (f>p\ (assuming (po > <!>*) 
as 



Ci = -t 



ypi 



)) = -!• 



(A5) 



Finally, Eq. (jA4p can be integrated again; the general 
solution can be written as 



C2 



bo - (j)i)fi{(j)i)d(j)i 

(A6) 



where C2 is an integration constant, and we introduced 
the auxiliary function 



m(</') = exp 



V 



exp 



3Af4j / 



\V{(b.) 1/(0) 



(A7) 



The value of C2 is determined through the boundary con- 
dition at = 0, . However, we note that /z(0) grows 
rapidly with 0; therefore the C2 term is negligible for 
away from the reheating point. The closed-form expres- 
sion for ^ is thus 



1 



D{(b)ti{cP) J^^ ^ 

-j^ ^min(0,0o) 



bo - 4>i)^J.{(j)i)d(j)i 



t>i exp 



D{ 



d(j>2 



(A8) 



We will now derive simplified forms of this expression in 
cases (j) < (j)o and 4>> 4>o- 

When < 00 (but for not too close to 0*), we can 
simplify Eq. (jASp if we note that the outer integrand is 
dominated by the neighborhood of 0i = where the 
exponent is close to 1. Then we can perform an asymp- 
totic estimate of the integral in Eq. (|A8p . The easiest 
method is to integrate by parts repeatedly, which yields 
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an asymptotic series: 



I'D, ' 



1 



D{<P) 
1 



D{4>i) 



exp 



exp 



D 



D ^ 



D{^) , D{4>) (D{^)\ 



Di^) vicb) v{(b) V «W 



(A9) 



In retaining only the first term of the asymptotic series, 
we neglect terms involving H/Mpi <?C 1 as well as terms 
proportional to the slow-roll parameters. 



For completeness, we give the result also for cj) > (f)o. 
In that case, the integral over in Eq. (|A8[) becomes 
(/)-independent, and we get 



Here we neglected the exponentially small terms of order 



exp 



r4> 


V 




1 




D 




m(0) 



< 1. 



Thus we find for (p < (pQ the required result, 



(AlO) 



1 



exp 



8 \V{(Po) V{<P) 



(A12) 



v{4>) 



Dm ' 



v{cP) 
(All) 
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